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We perform both analytical and numerical studies of the one-dimensional tight-binding Hamilto- 
nian with stochastic uncorrelated on-site energies and non-fluctuating long-range hopping integrals 
Jmn ~ J/\m — It was argued recently [A. Rodriguez at al, J. Phys. A: Math. Gen. 33, L161 
(2000)] that this model reveals a localization-delocalization transition with respect to the disorder 
magnitude provided 1 < fi < 3/2. The transition occurs at one of the band edges (the upper one 
for J > and the lower one for J < 0). The states at the other band edge are always localized, 
which hints on the existence of a single mobility edge. We analyze the mobility edge and show 
that, although the number of delocalized states tends to infinity, they form a set of null measure in 
the thermodynamic limit, i.e. the mobility edge tends to the band edge. The critical magnitude of 
disorder for the band edge states is computed versus the interaction exponent jj, by making use of 
the conjecture on the universality of the normalized participation number distribution at transition. 

PACS numbers: 71.30.-|-h; 72.15.Rn; 78.30.Ly; 36.20.Kd 



I. INTRODUCTION 

In 1958 Anderson formulated a simple tight-binding 
model with uncorrelated on-site (diagonal) disorder and 
predicted a localization-delocalization transition (LDT) 
in three dimensions (3D): a phase of extended states ap- 
pears at the band center in the thermodynamic limit if 
the disorder magnitude is smaller than a critical value, 
while all states are localized at larger magnitudes of dis- 
orderji The phase of delocalized states is separated from 
the two phases of localized states by two mobility edgesjS 
The concept of the mobility edge is of key importance for 
the low-temperature transport properties of disordered 
materials. 

Since the advent of the single-parameter scaling hy- 
pothesis, introduced by Abrahams et al.,^ the occurrence 
of a localization-delocalization transition (LDT) in disor- 
dered systems with time reversal symmetry was ruled out 
in one- (ID) and two dimensional (2D) geo metr ies at any 
disorder strength (for the overview see Refs. l4l5l lgl7l). The 
localized nature of the states in ID was pointed out even 
earlier by Mott and Twosei^ 

At the end of eighties and beginning of nineties 
it was realized however that correlations in disor- 
der may give rise to extended states in low dimen- 

Sions,9,l"i".12il3,14,15,16,17,18,19,20,21,22,23,24 r^^^g^ ^-^Ovt- 

range correlated on-site disorder was found to cause the 
appearance of extended states at special resonance ener- 
gies in ID. They form a set of null measure in the density 
of states in the thermodynamic limit jSiiStiiiiSiiiM imply- 
ing the absence of the mobility edges in those models. In 
spite of this fact, even the infinitesimal fraction of the 
extended states may have a strong impact on the trans- 
port properties of disordered materials. In particular, 
short-range correlations in disordered potential were put 



forward to explain unusual conducting properties of poly- 
mers, such as polyaniline and heavily doped polyacety- 
lene^^^ as well as semiconductor superlattices grown 
with random but correlated quantum well sequences^ 

Contrary to short-range correlations in the disorder 
distribution, long-range correlations were demonstrated 
to cause the LDT in ID systems, which is analogous 
to the standard Anderson LDT in 3D,i6,i7ii8,i9^20,2i 
this regard, a ID system with nearest-neighbor interac- 
tions and long-range correlated on-site disorder distribu- 
tion with a power-like spectrum S{k) ^ k^°' is critical 
with respect to the exponent a. More specifically, when 
the standard deviation of the energy distribution equals 
the nearest-neighbor hopping and a < 2 all states are 
localized, while for a > 2 a phase of extended states ap- 
pears at the center of the band giving rise to two mobility 
edges. The phase occupies a finite fraction of the density 
of states. A similar picture holds in 2Dm^ The authors 
of Ref. 18 proposed to use the long-range correlated dis- 
order and the appearance of a phase of extended states 
for designing microwave filters based on a single-mode 
waveguide. This type of disorder is also being studied in 
biophysics in connection with the large-distance charge 
transport in DNA sequences. ^^'^^ 

Another ID model which exhibits the LDT and a phase 
of extended states is an ensemble of power-law random 
banded matrices Hik cx Gik \i — where the matrix 

Gik runs over a Gaussian orthogonal ensemble2L2& (for 
an overview see Ref. '2^ . This model is critical with re- 
spect to the interaction exponent a: for a > 1 all states 
are localized, while all of them arc delocalized at a < 1, 
suggesting that a = 1 is the critical point in the model. 
Within the framework of this model, it was demonstrated 
rigorously that (i) the distribution function of the inverse 
participation ratio is scale invariant at transition and 
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(ii) the relative fluctuation of the inverse participation 
ratio (the ratio of the standard deviation to the mean) is 
of the order of unity at the critical point i^S*^ This find- 
ing confirmed the conjecture (which was put forward ear- 
lieri^^i^), that distributions of relevant physical magni- 
tudes are universal at criticahty (see also Refs. ls^ls^IlM) . 
The invariance can therefore be used to monitor the crit- 
ical point. •^'''■^^ 

Recently, several reports addressed to the unusual lo- 
calization properties of ID and 2D tight-binding mod- 
els with uncorrelated diagonal disorder and non-random 
long-range coupling between sites m and n, which falls 
according to a power law J/|m — n|''(2L22i22ii2iiiii^ (see 
also Refs. I43ll4# . More specifically, the states at one of 
the band edges (the upper one for J > and the lower 
one for J < 0) undergo the LDT with respect to the disor- 
der strength A if the interaction exponent /i ranges within 
the interval 1 < < 3(i/2, d being the dimensionality4i 
In what follows we set J > 0, so that extended states can 
appear at the upper band-edge. The states at the other 
band-edge are strongly localized, no matter how small the 
disorder magnitude is, thus, suggesting the existence of a 
single mobility edge.''*' At ^ > 3d/ 2 all states were found 
to be localized. The character of localization, however, 
turned out to be governed by the interaction exponent 
^. At ^ = 3d/2, the upper band-edge states are weakly 
localized, similar to those at the center of the band in the 
standard 2D Anderson model)^ while at fi > 3d/2 they 
are strongly localized. 

In the present paper we deal with the analysis, both 
numerical and analytical, of the localization properties 
of the latter class of ID Hamiltonians. We calculate the 
phase diagram of the transition for the upper band-edge 
states, i.e., the critical disorder magnitude Ac versus the 
interaction exponent /i, and show that Ac vanishes as 
(3/2 — /i)^/'^ when fi 3/2, while it diverges as 
when fi 1. Applying the finite size scaling analysis we 
study the problem of the mobility edge and show that 
the mobility edge approaches to the upper band edge in 
the thermodynamic limit. In other words, the fraction 
of the delocalized states forms a set of null measure, al- 
though their number tends to infinity on increasing the 
system size N as TV^^/^"''^/^^"'') (note that the depen- 
dence is sub-linear). This numerical finding is supported 
by a simple qualitative arguments based on the compar- 
ison of size scaling of two magnitudes: the bare level 
spacing at the band-edge and effective disorder seen by 
the quasiparticle. 

The outline of the paper is as follows. In the next Sec- 
tion, we describe the model and briefiy overview quali- 
tative arguments which brought us to the conjecture on 
the existence of the LDT within the model. In Sec. IIIII 
we present the phase diagram of the LDT for the up- 
per band-edge states, which is calculated on the basis of 
the statistics of the participation number. The mobility 
edge and fraction of the delocalized states are addressed 
in Sec. II VI We conclude the paper in Section Ivl 



II. MODEL AND QUALITATIVE REASONING 

We consider a tight-binding Hamiltonian on a ID reg- 
ular lattice with N sites 



H = ^e„|n)(n| + ^ J,nn\m){n\ 



(1) 



Here \n) is the ket vector of a state with on-site energy 
£„. These energies are taken at random and uncorre- 
lated for different sites and distributed uniformly around 
zero within the interval [—A/2, A/2], having therefore 
zero mean, (e„) = 0, and standard deviation (e^)'^^^ = 
A/ (the angle brackets (• • • ) denote average over dis- 
order realizations). The hopping integrals Jmn do not 
fluctuate and are set in the form Jmn = J/\m — n\^, with 
J > and Jnn = 0. 

First, we address the disorder- free system (A — 0), 
taking periodic boundary conditions for the sake of sim- 
plicity. Then the eigenstates of the Hamiltonian |^ are 
plane waves with quasi- momenta K = 27rfc /N within the 
flrst Brillouin zone k E [—N/2, N/2). The corresponding 
eigenenergies are given by 



n>0 ' ' 



(2) 



where the summation runs over all N sites of the lat- 
tice, and is assumed to be larger than unity to get a 
bounded energy spectrum. The complete account for all 
terms in the sum is important in the neighborhood 
of the upper band-edge, where the long-range hopping 
terms affect the dispersion vastly. At the upper band- 
edge {K 0) the dispersion relation ^ is as follows^ 

El°P{K) = ^top(M) - JAtop(/i)|i^r"' - JBt,p{ii)K^ , 

(3) 



when M 7^ 3. Here, Eto■p{^i) = '^J^,'^ 



2JC(^) is 



the upper band-edge energy in the thermodynamic limit, 
^top (m) 8.nd i?top (m) are dimensionless positive constants 
on the order of unity (for brevity we do not provide ex- 
plicit expressions for them), and C(/j,) is the Riemann 
^-function. For small K, the sub-quadratic term in the 
right-hand side of Eq. ||2J) dominates over the quadratic 
one if /i < 3 and vice versa. The range fi < 3 will be 
of our primary interest. We will focus later on the size 
scaling of the energy spacing at the upper band edge at 
/i < 3, which is 



SE'°P oc iV^-^ 



(4) 



At the lower band edge {\K\ tt), the energy spectrum 
is parabolic: 



E^°^K) = E^M + JB,,M (^ - KY 



(5) 



where i?bot(M) = 2 JYl^i^^)^''^^^ the lower band en- 
ergy, which depends weakly on /x, and -Bbot(M) is a di- 
mensionless constant on the order of unity. Respectively, 
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the energy spacing at the bottom of the band scales as 



1000 



5E^°^ cx N- 



(6) 



On introducing the disorder the eigenstates of the reg- 
ular system couple to each other, which can result in 
their localization. The typical fluctuation of the coupling 
matrix, namely the first term on the right-hand side of 
Eq. 1^ in the i^-space basis, is^ 



(7) 



^/VlN 



Now, we compare the size dependence of a to that of the 
eigenenergy spacing at the top of the band SEj^^^ (defined 
by Eq. |0J). As the system size increases the typical value 
of the eigenstate coupling a decreases faster than the en- 
ergy spacing SEj^P at ^ < 3/2. Consider a system of size 
N and let a be of perturbative magnitude (i.e. a <C SE), 
then the upper band-edge states are only weakly per- 
turbed by disorder and remain extended over the whole 
system. Upon increasing the system size, the inequality 
a SE gets even stronger, so that the perturbation be- 
comes weaker and therefore the upper states will remain 
delocalized at TV ^ oo. On the other hand, large disorder 
(say, larger than the bare bandwidth) would certainly lo- 
calize all the states. These arguments indicate that there 
can exist extended states at the upper band-edge at finite 
disorder, provided that 1 < /i < 3/2. This conjecture was 
confirmed both theoretically, by means of the renormal- 
ization group approach combined with a super-symmetric 
method for disorder averaging^, and numericallyi^£i^2il£ 
The value ii — 3/2 represents the marginal case in which 
the upper band-edge states are localized weaklyi^L^Siii 
In what follows, we focus therefore on the interaction 
exponent /i ranging within the interval 1 < < 3/2. 

At the bottom of the band {\K\ — > tt) the level spacing 
diminishes as 7V~^, that is faster than the effective magni- 
tude of disorder a. Thus, even if cr ^ SE for some lattice 
size and the states are extended over the whole system, 
the inequality will be reversed for larger TV, which will 
finally result in strong coupling of the states and their 
eventual localization. The same conclusion holds for the 
entire energy spectrum if fi > 3/2. 

The above picture implies the existence of a single mo- 
bility edge separating the phases of localized and delocal- 
ized states. We address this question in detail in Sec. lIVI 



III. PHASE DIAGRAM OF THE TRANSITION 

In this Section we calculate the dependence of critical 
disorder magnitude Ac on the interaction exponent /i G 
(1,3/2]. To detect the transition we analyze the wave 
function statistics. More specifically, we calculate size 
and disorder dependencies of the relative fluctuation of 




FIG. 1: Critical magnitude of disorder Ac as a function of the 
interaction exponent fi. Full circles represent the numerical 
data, while solid lines are the best fits: (a) — > 1 and (6) 
fj, -» 3/2. 



the participation number (PN) defined as 

-1 



Pk= (EI^^"!") 



(8) 



where ipkn denotes the n-th component of the normalized 
k-th eigenstate of the Hamiltonian JQ. Here, the state 
index k ranges from I to iV. By definition, we adscribe 
band-edge states k — 1 and k — N to the uppermost and 
lowermost eigenstates, respectively. It was demonstrated 
recently'^^'^^ that within the considered model, the ratio 
of the standard deviation of the PN (SDPN) to the mean 
of the PN (MPN) is an invariant parameter at transition 
and therefore can be used to detect the critical point. We 
note that wave function statistics turned out to be more 
efficient for this purpose than level statistics (the latter 
also represents a method to monitor the transition^''), 
because wave function statistics appeared to be less af- 
fected by strong finite size effects within the model (for 
details see the discussion in Refs. I37ll33). 
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As the LDT occurs at the top of the band within the 
considered model, we calculated disorder and size depen- 
dencies of the ratio SDPN/MPN for the uppermost state. 
Open chains were used in all calculations. Following the 
procedure developed in Ref. .38. we obtained the critical 
magnitude of disorder Ac at which uppermost eigenstates 
undergo the LDT for fi e (1,3/2]. The results of the sim- 
ulations are shown in Fig.^by the full circles, while solid 
lines represent best fits close to the limiting points fJ- — 1 
(Fig.Ch) and fi = 3/2 (Fig.^). We found that 

Ac « 10.1J(/i- 1)^1 , /^^l, (9a) 

Ac w32.6J(3/2-^)°-^V /i^3/2, (9b) 

First of all, we notice that, according to Eq. H9a|l . the 
critical magnitude of disorder diverges as Ac oc (/i — 1)^^ 
when /X ^ 1. The explanation of the divergence relies 
in the fact that Ac oc EJ^p for Ac > 1^ and the up- 
per band-edge energy E^°]^-^ diverges when iV ^ oo as 
{fi — 1)~^. Contrary to that, the critical magnitude of 
disorder vanishes as Ac oc (3/2 — /x)"-^^ when /i 3/2, 
indicating that in the marginal case (/i = 3/2) the up- 
permost state is localized (Ac = 0), in full agreement 
with the results obtained by the renormalization group 
approach combined with a super-symmetric method for 
disorder averaging 

IV. MOBILITY EDGE AND FRACTION OF 
THE DELOCALIZED STATES 

Having discussed the dependence of the critical mag- 
nitude of disorder Ac on the interaction exponent fi for 
the uppermost state, we now focus on the mobility edge 
and the fraction of extended states in the thermodynamic 
limit. This question has not been addressed in previous 
studieS)^£i22*^2iiMiiiS except for a comment on the exis- 
tence of a single mobility edge^Si^ (see also the discussion 
in Secini). 

To work out this problem, we numerically diagonalized 
the Hamiltonian ^ and calculated the normalized MPN, 
{P)/N, as a function of energy for different system sizes 
N. The results of these simulations are shown in Fig. 13 
for a particular set of interaction exponent and magni- 
tude of disorder {fi — 5/4 and A — 5J). First, one can 
see that the upper band-edge energy increases with N as 
was mentioned in Sec. lIIII fsee also Refs. i37«38i') . Note also 
the noticeable size dependence of the normalized MPN 
of the uppermost state. This behavior reflects the finite 
size effects already mentioned above (see the preceding 
section). The boundaries result in a positive correction of 
the order of N~'^ to the bulk value of the uppermost wave 
function (~ l/^/W). Consequently, the normalized MPN 
depends on the system size as (1 — cN^^f^), where c is a 
constant. The second, and the most important, observa- 
tion is that the normalized MPN increases monotonically 
on approaching the upper band edge, for all considered 
values of N. 




E/J 



FIG. 2: The normalized mean participation number {Pk)/N 
as a function of energy for different system sizes calcu- 
lated for the interaction exponent fi — b/i and the disorder 
strength A = 5J. 




k/N 

FIG. 3: The mean participation number as a function of the 
normalized state index k for different system sizes A'' calcu- 
lated for the interaction exponent /i = 5/4 and the disorder 
strength A = 5J. The inset shows an enlarged view close to 
the upper band edge. 



In order to avoid the size dependence of the upper 
band-edge, we use hereafter the state index k rather than 
the energy Ef^{k). In Fig. 13 we plotted the MPN as a 
function of the normalized state index k/N. The fig- 
ure demonstrates that the MPN is independent of the 
system size in a wide range of the normalized state index 
k/N. The perfect collapse of the curves within this range 
clearly indicates the localized nature of these eigenstates. 
However, at the top of the band the MPN increases lin- 
early with the system size and the collapse is absent (see 
the blow-up in the inset of Fig.|2|). This result suggests 
that not only the uppermost eigenstate, but a number 
of them are delocalized, in agreement with the previous 
claim raised in Ref. i41| 
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We now apply a finite size analysis of our numerical 
results to further confirm the latter statement and to 
obtain the size scaling of the number of extended eigen- 
states. Figure 0] shows that the MPN for different sizes 
collapses onto a single curve close to the top of the band 
after introducing the rescaled index k/N^/"^ for /i = 5/4. 
The collapse holds up to a finite value of the rescaled in- 
dex, which indicates that the number of extended states 
scales as oc N^^^ when ^ — 5/4. 




FIG. 4: Normalized mean participation number as a function 
of k/N^^^ for different system sizes calculated for fi = 5/4 and 
A = 5J. 

To provide support to these numerical results, we now 
develop a simple (pcrturbative) analytical approach that 
allows us to understand the origin of the obtained results. 
To this end, we analyze the energy spacing close to the 
upper band-edge in the disorder-free system (A = 0) 
more accurately. The spacing can be obtained from 
Eq. ®: 



SEl°Pik) = JC, 



1 



(10) 



The constant absorbs all unessential numerical fac- 
tors. Note the absence of the scaling of the spacing 
(5i?*°P(fc), i.e., the spacing depends on k not only on k/N. 
This is crucial for the peculiar features of the LDT within 
the model. Indeed, consider a chain of a particular size 
N. Assume that the disorder is pcrturbative for the up- 
permost state (fc = 1), i.e., a ^ SEj^P{l), so that states 
at the top are extended over the whole sample. Let us 
now find the mobility edge for this finite lattice, defining 
it by the equality 



SEl^^ik, 



(11) 



The number km, which is A^-dependent, divides all eigen- 
states into two sets: those with k < km are delocalized 
in the above sense, while eigenstates with k > km are lo- 
calized in the usual sense. Thus, km provides us with the 



number of extended states for a particular system size 
N. Determined by the equation it reads 



km, ^ ( 



1/(2-m) 



7V(3/2-m)/(2-m) ^ (12) 



where all unessential constants are absorbed into D^. 
Applying Eq. p2(l to the particular case /i — 5/4 we 
recover the behavior that we have found numerically, 
namely km oc N^^"^. 

The relationship (|12() provides us also with the depen- 
dence of km on the disorder magnitude A, which can also 
be compared to numerical calculations. We performed 
such comparison for the interaction exponent /i = 5/4. 
In this particular case, the exponent 1/(2 — /i) in the A- 
dependence of km is equal to 4/3. Figure El shows the 
normalized mean participation number as a function of 
fcA^/'^ calculated for different disorder strength A and a 
given system size A^ — 4000. The collapse of all curves 
onto a single one in the vicinity of the upper band-edge 
supports the validity of Eq. ((T^ . 
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FIG. 5: Normalized mean participation number as a function 
of kA^^^ for different disorder magnitudes with = 5/4 and 
N = 2000. 

We have also tested Eq. (|12(l against numerical simu- 
lations performed for other values of /j, G (1,3/2) and al- 
ways found an excellent agreement of the numerical data 
with the analytic formula p2|l . It follows from the for- 
mula that the number of extended states increases sub- 
linearly with the system size, namely as iV(3/2-^)/(2-^)^ 
which means that the fraction of these states vanishes as 
]\[-i/i^~^tJ-) when A^ — > oo. From this we conclude that 
the mobility edge approaches to the upper band edge in 
the thermodynamic limit. 



V. SUMMARY AND CONCLUDING REMARKS 

We studied both analytically and numerically the lo- 
calization properties of the ID tight-binding model with 
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diagonal disorder and non-random long-range inter-site 
interactions, Jmn — J/\'m — n\^ where J > 0. The 
model can be critical at the upper band-edge provided 
1 < < 3/2. 

We calculated the phase diagram of the transition (the 
dependence of the critical magnitude of disorder Ac on 
the interaction exponent ^) by studying the participation 
number statistics. The critical magnitude of disorder was 
detected using the size invariance of the ratio of the stan- 
dard deviation to the mean. We found that diverges 
as {jjL—l)^^ when /i ^ 1, which originates from the simi- 
lar divergence of the upper band-edge energy. If /i 3/2, 
the critical disorder magnitude vanishes as (3/2 — fi)"^^^, 
indicating that all states are localized at n = 3/2, no 
matter how small the disorder is. 

It is shown, both analytically and by means of the fi- 
nite size scaling analysis, that the number of extended 
states at the upper band-edge increases sub-linearly with 
the system size N, namely as (x iV^^/^"'^^/^^"^)), there- 
fore forming a set of null measure in the density of states 
in the thermodynamic limit. This suggests that the mo- 
bility edge is only a meaningful concept for a finite size 
system; it approaches to the upper band-edge in the ther- 
modynamic limit. Although extended states form an in- 
finitesimal fraction of the whole density of states, they 
can provide a strong impact on the transport properties. 



similarly to what happens in systems with correlated ran- 
domness J^iSii^ 

To conclude, we note that our findings, apart from be- 
ing interesting from the theoretical point of view, are rele- 
vant for real physical systems. Thus, some organic mate- 
rials with planar geometry, in which optically allowed ex- 
citations are dipolar Frenkel excitons ;'^^i'^^''^7i'^?i4?iWi^ii^'?^ 
represent an example; the value of the interaction expo- 
nent for the dipole-dipole inter-site interaction is ii — 3, 
which resembles the weak localization regime in the ID 
long-range model when /i = 3/2.^^ Dipole-exchange spin 
waves in ferromagnetic films provide yet another example 
where these results are relevant 
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